Optimal preparation of graph states 
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We show how to prepare any graph state of up to 12 qubits with: (a) the minimum number 
of controlled-Z gates, and (b) the minimum preparation depth. We assume only one-qubit and 
controlled-Z gates. The method exploits the fact that any graph state belongs to an equivalence 
class under local Clifford operations. We extend up to 12 qubits the classification of graph states 
according to their entanglement properties, and identify each class using only a reduced set of 
invariants. For any state, we provide a circuit with both properties (a) and (b), if it does exist, 
or, if it does not, one circuit with property (a) and one with property (b), including the explicit 
one-qubit gates needed. 

PACS numbers: 03.65. Ud, 03.65. Ta, 03.67.Mn, 42.50.Xa 



I. INTRODUCTION 

Building a quantum computer entails isolating atomic- 
scale systems, except when a controlled interaction (e.g., 
a logic gate) is applied. Isolation must be kept up until a 
subsequent controlled interaction is applied. Any unde- 
sirable coupling with the outside disrupts the quantum 
state of the systems and ruins the computation. Due 
to the enormous difficulties which keeping atomic-scale 
systems isolated and controlled entails, a main limiting 
factor preventing the development of quantum comput- 
ing is the amount of time during which the systems must 
be kept isolated and controlled to achieve the computa- 
tion. Therefore, the task of reducing this amount of time 
is an essential factor in the process of achieving the goal 
of building a quantum computer. 

Another limiting factor for quantum computation is 
the number of entangling gates needed. While one-qubit 
gates can be built with fidelities higher than 99%, two- 
qubit entangling gates hardly reach 93%, and this be- 
comes worse for three-qubit gates, etc. Therefore, since 
one- and two-qubit gates are enough for universal quan- 
tum computation, it is reasonable to focus on circuits 
with only one- and two-qubit gates, and having the least 
possible number of two-qubit gates. 

Here we address the problem of preparating an im- 
portant class of quantum states with circuits requiring 
(a) the minimum number of two-qubit entangling gates, 
and (b) the minimum preparation depth (i.e., requiring 
a minimum number of time steps). We assume that we 
can implement arbitrary one-qubit gates and one specific 
two-qubit entangling gate, the controlled-^ gate. The 
result can be easily extended to any other specific two- 
qubit gate. 
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A. Graph states 

Graph states [l|, Q are a type of n-qubit pure state 
with multiple applications in quantum information pro- 
cessing. Two important examples are: in quantum error- 
correction, the stabilizer codes which protect quantum 
systems from errors [H-IH and, in measurement-based (or 
one-way) quantum computation M , the initial states con- 
sumed during the computation [7[. 

The definition of a graph state already provides a 
blueprint for its preparation: an n-qubit graph state \G) 
is a pure state associated with a graph G = (V, E) con- 
sisting of a set of n vertices V = {0, . . . , n — 1} and a 
set of edges E connecting pairs of vertices, E C V x V 
[HQ. Each vertex represents a qubit. An edge G E 
represents an Ising-type interaction between qubits i 
and j. To prepare |G), first prepare each qubit in the 
state |+) = 4g(|0) + |1))> i- e -j the initial state will be 
IV'o) = 0iev |+)»- Then, for each edge G E con- 

necting two qubits i and j, apply a controlled- Z gate 
between these qubits, i.e., the unitary transformation 
Cz = |00)(00| + |01)(01| + |10}(10| - |11)(11|. 



B. Preparation using only controlled- Z gates 



c have the state |^o) = 
to prepare the eight-qubit 



Let us suppose that 1 

®iLo an d wc want 

graph state \G) whose graph G is in Fig. [TJ using only 
controlled- gates. One possible way is to follow the 
preparation procedure suggested by G, taking into ac- 
count the possible restrictions in performing two or more 
controllcd-Z gates simultaneously, and optimally dis- 
tributing the controlled- Z gates to minimize the number 
of steps. The preparation depth of \G) is the minimum 
number of time steps required to prepare \G) Q. 

The state \yj ) corresponds to a graph with n isolated 
vertices. The total number of edges in G gives a trivial 
upper bound to the preparation depth of \G), since each 



FIG. 1. Graph corresponding to the 8-qubit graph state we 
want to prepare. Edges in the same color represent controlled- 
Z gates that can be performed in the same time step. 



controlled- can be applied in a time step. To find the 
minimum number of time steps, we have to explore the 
possibility of applying two or more controlled- Z gates in 
a single time step. 

Given a vertex i in G, the neighborhood of i, 
is the set of vertices connected to i. Now let us sup- 
pose that there is more than one element in Af(i), i.e., 
\N(i)\ > 1, and let j, k £ Af(i) be two of the neighbors. 
Then, to prepare the corresponding graph state \G) we 
must apply a controlled- .Z gate to entangle qubits i and 
j, and another one to entangle qubits i and k. Applying 
a controlled- Z gate between qubits i and j in a certain 
time step implies that both qubits are busy during this 
entangling interaction. If we focus on qubit i, we have 
to wait for a further time step to have qubit i free be- 
fore applying another controlled- Z gate to entangle i and 
k. Nevertheless, vertex k could be connected to another 
vertex I ^ in G. If such is the case, we could in 

principle take advantage of the same time step we are 
using to entangle qubits i and j in order to do the same 
with qubits k and I, because both controlled- Z gates can 
be performed in parallel. 

Remarkably, the problem of determining the minimum 
number of time steps with the restrictions that we have 
mentioned is related to an old problem in graph theory: 
given an edge in G, let us use a certain color to 

mark and those other edges of G corresponding to 
controlled- gates that can be performed at the same 
time step than that of and use different colors for 

those edges related to controlled- Z gates that do not ful- 
fill this condition. Since two controlled- Z gates can be 
performed at the same time step if and only if the four 
qubits involved do not coincide, then every edge incident 
to the same vertex must have a different color. In graph 
theory this color configuration is called a proper edge col- 
oring or, to put it more concisely, the graph is said to be 
edge-colored. Hence, the preparation depth problem is 
equivalent to determining the minimum number of col- 
ors required to get a proper edge coloring of G. This 
problem is known as the determination of the chromatic 
index or edge chromatic number ~x!(G). 

Let us denote by A(G) the maximum degree of G (i.e., 
the maximum number of edges incident to the same ver- 
tex). Vizing's theorem Q states that G can be edge- 



FIG. 2. A circuit with minimum preparation depth for the 
graph state corresponding to Fig. [T] using only controlled- Z 
gates. In the circuit, a qubit (vertex of the graph) is repre- 
sented by a horizontal wire, and a controlled-/? gate (edge in 
the graph) by a vertical segment with diamond-shaped ends 
connecting the qubits involved in the operation. Time steps 
are separated by vertical dashed lines. 



colored in either A(G) or A(G) + 1 colors, and not fewer 
than that. Therefore, x'(G) is either A(G) or A(G) + 1. 
A proof can be found in [lfj . The important point is that 
x'(G) gives the preparation depth of |G) (l.flil]. 

Graphs requiring A(G) colors are called class-1 graphs, 
and those requiring A(G) + 1 colors are called class- 
2 graphs. For instance, the four-qubit fully connected 
graph state is represented by a graph of four vertices, six 
edges, and maximum degree equal to 3: it is a class-1 
graph, so that its preparation depth coincides with its 
maximum degree: 3. On the other hand, the three-qubit 
fully connected graph state is represented by a graph 
of three vertices with three edges and maximum degree 
equal to 2: it is the smallest class-2 graph and, as a con- 
sequence, its preparation depth is also 3. 

The graph in Fig. [T] is a class-1 graph: it can be 
edge-colored with A(G) = 7 colors (this number corre- 
sponds to the degree of vertex 3). Hence, if we use only 
controlled- Z gates, the minimum preparation depth of 
the corresponding graph state is 7. For instance, one of 
the optimal distributions of the 13 controlled- Z gates is 
given in the circuit of Fig. [5J which has seven time steps. 

However, if we are allowed to use one-qubit gates, in 
most cases it is possible to get the desired graph state 
with a lower number of two-qubit gates and less prepa- 
ration depth. For that purpose, one has to take into 
account the degree of entanglement of the state we want 
to prepare. 



C. Optimum preparation 

Two graph states have the same degree of entangle- 
ment if and only if they are equivalent under local uni- 
tary (LU) operations [12| , so that they belong to the same 
LU-equivalence class. Moreover, previous results suggest 
that, for graph states of up to 26 qubits [l3| . the no- 
tions of LU-equivalence and LC-cqui valence (equivalence 
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under local Clifford operations) coincide and, therefore, 
entanglement classes are in fact local Clifford equivalence 
classes (LC classes). This implies a remarkable simplifi- 
cation, since it is possible to carry out a graphical de- 
scription of the action of local Clifford transformations 
on graph states [l4| : the successive application of a sim- 
ple graph transformation rule, the so called local comple- 
mentation, on a certain graph G and on those that are 
obtained from G, allows us to generate the whole LC class 
of |G). The entire set of graphs connected to a given G by 
a sequence of local complementations is usually referred 
to as the orbit of G (i.e., the LC class of \G)). 

Any of the graph states belonging to a given LC class 
could be used as a representative for that orbit, but there 
is a practical advantage in taking one requiring: (a) the 
minimum number of controlled- gates, or (b) the min- 
imum preparation depth of the class. If we can identify 
which LC class a given \G) belongs to, then we can pre- 
pare \G) by preparing instead the LC-equivalent state 
G") requiring (a) and (b), and then transform \G') into 
\G) by means of single-qubit Clifford operations corre- 
sponding to a sequence of local complementations car- 
ried out on G' . This last transformation from \G') to 
\G) requires only one additional time step, since the lo- 
cal complementation at vertex a is equivalent [l], Q to 
the unitary operation 

U:(G)=exp(-i^) J] exp^W), (1) 

where cri°^ is the Pauli matrix a x acting on qubit a, 
is the Pauli matrix a z acting on qubit b, and Af(a) is the 
neighborhood of a. Operations on different qubits com- 
mute, and therefore one can group together a sequence of 
local complementations UJ L (G) into n one-qubit gates Ri 
(with i = 0, . . . , n — 1), which can be jointly performed 
in a single step. 

The preparation procedure of a graph state \G) 
through the optimum LC-representative provides an ad- 
vantage in time steps when compared to the standard 
graph-based controlled- Z procedure when 

X , (G)- X '(G')>1, (2) 

since the preparation depth for the standard procedure is 
x'(G), while the preparation depth through the optimum 
LC-representative is x'(G') + 1, that is, the sum of the 
preparation depth of the optimum LC-representative plus 
a unit of time corresponding to the one-qubit gates. This 
means that the preparation depth through the optimum 
LC-representative provides an advantage for most graph 
states. For instance, as we will see, for the graph state 
of Fig. [TJ the optimum preparation circuit requires only 
seven controlled- Z gates and three time steps: saving six 
controlled-.!? gates and requiring four time steps fewer 
than in the standard procedure. 

However, the preparation through the optimum LC- 
representative requires us to identify which LC class \G) 



belongs to. To be practical, this must require us to iden- 
tify the simplest signature of the class. 

In Sec. mi we classify all LC classes for graph states 
up to n — 12 qubits. This classification is based on a re- 
duced set of invariants which allows us to identify which 
class a given state belongs to. In Sec. IIII1 we provide a 
representative of the class with both properties (a) and 
(b), if it exists, or, if it does not, one with property (a) 
and one with property (b). All these results, which oc- 
cupy several hundreds of megabytes, are organized in two 
tables, one for n < 12 and one for n = 12, and presented 
as supplementary material [15| . Finally, in Sec. IIV[ we 
explain how to obtain the one-qubit gates needed, and 
provide as supplementary material a computer program 
[lq to, given the graph G corresponding to the state we 
want to prepare, generate a sequence of local complemen- 
tations which connect G to the corresponding optimum 
graph(s). In Sec. El the whole method is illustrated with 
an example. 



II. CLASSIFICATION OF GRAPH STATES IN 
TERMS OF ENTANGLEMENT 

The classification of the entanglement of graph states 
has been achieved up to n = 7 [l], , and has recently 
been extended to n = 8 [TtJ ■ There are 45 LC classes for 
graph states up to seven qubits, and 101 LC classes for 
cight-qubit graph states, which are ordered according to 
certain criteria based on several entanglement measures, 
which arc invariant under LC transformations. 

Here we extend the classification up to n = 12. The 
number of orbits for n < 8 qubits is 146. For n = 9, 
there are 440 orbits; for n = 10, there are 3 132 orbits; 
for n — 11, there are 40 457; and, for n = 12, there are 
1 274 068. All the information about each LC class is 
presented as supplementary material fl5| . 

For each LC class we give the number of nonisomor- 
phic graphs in that class (size of the orbit), \LC\. Then, 
the orbits are classified according to the number of ver- 
tices (qubits), |V|; the minimum number of edges of a 
graph belonging to the class (controlled- Z gates needed 
for the preparation), \E\; the Schmidt measure, E$ (up- 
per and lower bounds are given where the exact value is 
unknown); and for n/2 > i > 2, the rank indexes RI{ for 
bipartite splits with i vertices in the smallest partition. 
The classifying labels |V|, \E\, E$, and Rli are applied 
in this order. Additionally, we include the information 
regarding the existence (or not) of a two-colorable graph 
belonging to the class (a piece of data which is useful, in 
some cases, to calculate lower and upper bounds for the 
Schmidt measure). 

However, these numbers, which were (almost) enough 
to identify every class if n < 8, are not enough to identify 
every class if n > 8. In other words, the set of entangle- 
ment measures for n-qubit graph states used in [H l2TTl7T| 
failed to distinguish between inequivalent classes under 
local Clifford operations if n > 8: different LC classes 
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had coincident values for those entanglement measures. 
In fact, the number of problematic LC classes increases 
with n. Therefore, using these invariants for deciding 
which entanglement class a given state belongs to is un- 
reliable. A finite set of invariants that characterizes all 
classes has been proposed in [l8j]. However, already for 
n = 7, this set has more than 2 x 10 36 invariants which 
are not explicitly calculated anywhere, and hence this set 
is not useful for classifying a given graph state. 

Nevertheless, a compact set of invariants related to 
those proposed in [l8| is enough to distinguish among 
all inequivalent LC classes with n < 8 qubits: the 
4 two-index invariants called cardinalities-multiplicities 
fl9l ] . There is a straightforward procedure for calculating 
these four invariants using the information contained in 
the graphs. 

We have analyzed the utility and limitations of the 
cardinalities-multiplicities (C-M hereafter) as LC-class 
discriminants for graph states up to 12 qubits, a ques- 
tion that was left as an open problem in [l9j . where it 
was conjectured that four of these C-M invariants would 
be enough to label and discriminate all the LC classes. 
Our results show that for graph states of n > 9 qubits the 
C-M invariants fail to distinguish between inequivalent 
LC classes: the smallest counterexample of the conjec- 
ture corresponds to a pair of nine-qubit orbits that have 
exactly the same entire list of C-M invariants. These 
are the only problematic orbits for n = 9 qubits. As an 
alternative for discriminating between them, we have cal- 
culated the whole list of Van den Nest-Dehaene-De Moor 
(VDD) invariants of type r = 1 [HI for these two orbits, 
and once again these invariants coincide. In order to 
determine the number of C-M and VDD "problematic" 
(undistinguishablc) orbits, we have extended our calcula- 
tions up to n = 12 qubits. The ratio Pf(n) of the number 
of graphs belonging to problematic orbits and the over- 
all number of graphs for each n, gives the probability 
that a randomly chosen graph state falls in one of the 
problematic orbits. The values of Pf(n) are, fortunately, 
quite low (see Table [I]). Therefore, the first step of the 
procedure of preparation, i.e., the identification of the 
orbit, resorts to C-M invariants (and, sometimes, type 
r = 1 VDD invariants), and works in most cases. For 
those rare states whose orbit identification through C-M 
or VDD (r = 1) invariants is not univocal, one would re- 
sort to VDD invariants of higher order r [l^]. However, 
the computational effort of this task makes this proce- 
dure less efficient than simply generating the whole LC 
orbit of the graph. 

In the supplementary material [HI, those LC classes 
which have the same set of quantities {\V\, \E\, E$, Rli} 
are ordered according to the C-M invariants, going from 
the smallest cardinality (zero) to the biggest one, and in- 
creasing the value of the associated multiplicity of each 
cardinality. For those orbits with the same set of quan- 
tities {\V\, \E\, Es, Rli, C-M invariants}, orbits are or- 
dered by the increasing size of the orbit, \LC\. We do not 
apply any subsequent classifying criteria, in case there 



TABLE I. Orbits for which C-M and VDD invariants are not 
good LC discriminants. 

n No. of orbits No. of problematic orbits Pf{n) 



< 



146 

9 440 

10 3 132 

11 40 457 

12 1 274 068 





2 
8 

78 
472 





0.0012218 
0.0006996 
0.0011929 
0.0000949 



were undistinguishablc LC classes left. 

C-M invariants are given as an ordered list of mul- 
tiplicities Mi, for i = 0,...,x. The value Mj is the 
multiplicity of the cardinality i. We do not list all 2 n 
possible multiplicities, but only the multiplicities of car- 
dinalities 0, . . . ,x, where x is the smallest number such 
that all "non-problematic" orbits are distinguished. This 
may not be the smallest possible set of C-M invariants. 
For instance, only the cardinalities {0, 1, 3, 4} are needed 
for n < 8, but, for simplicity's sake, we have included the 
continuous list {0, 1, 2, 3, 4}. 



III. OPTIMUM REPRESENTATIVE 

We have determined the optimum representatives for 
all the orbits up to 12 qubits, initially defined as one 
of those with the minimum number of edges in the or- 
bit and, among them, one of those with the minimum 
chromatic index. In the supplementary material (l5j . we 
have included a column labeled min(|i?|, x' , #) for each 
LC class, where \E\ is the minimum number of edges 
in the class, x' is the minimum chromatic index of the 
graphs with \E\ edges, and # is the number of noniso- 
morphic graphs with \E\ edges and chromatic index x' ■ 
The value of x' coincides with the preparation depth of 
those representative states and indicates how much more 
efficient the preparation procedure we are proposing is 
than the standard one [it is x'{G') in Eq. ©]. 

While carrying out the LC classification, an interest- 
ing observation arose: the optimum representative of a 
certain orbit was determined by the application of two 
filters to the orbit in a certain order. First, we looked 
for the graph with the minimum number of edges (which 
implies a minimum number of two-qubit entangling op- 
erations), \E\, and second the minimum chromatic index 
x' fixed that l-El (which means a minimum preparation 
depth given those graphs with \E\ edges). However, if we 
commute the order of the filters, for n < 7 qubits, the re- 
sult of the permutation of the filters gives the same graph, 
but this is not the case for n > 7. We get the simplest 
example of non coincidence between the final optimum 
representatives for n = 8, where there is only one orbit 
(LC class number 136) in which the permutation of the 
filters does not produce the same final graph. There are 
two graphs in this orbit with \E\ = 11 and x' = 4, and 
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one graph with x' = 3 and \E\ = 12 (see Fig. [3]). For 
each n < 12, we calculated the number of orbits with dif- 
ferent final representatives when we applied the filters in 
different order: the ratio of the number of "exceptional" 
orbits and the entire number of orbits increases with n, 
for n > 9 (Table [HJ. 



TABLE II. Orbits for which a different order of the filters re- 
lating the minimum number of edges and minimum chromatic 
index produces a non-coincident representative graph. 



# of orbits 



# of exceptional orbits 



< 7 
8 
9 
10 
11 
12 



45 
101 
440 
3 132 
40 457 
1 274 068 




1 

3 

65 
2 587 
136 518 



Moreover, we have included in our tables (see sup- 
plementary material [HI), beside the column labeled 
min(|i?|, x', #), another column labeled min(x', \E\, #), 
with a 3-tuple of values (%', \E\,#), where x' is now the 
minimum chromatic index in the class, \E\ is the mini- 
mum number of edges of the graphs with chromatic index 
x', and # is the number of nonisomorphic graphs with 
chromatic index x' and \E\ edges. Checking the coinci- 
dence of x' and \E\ in both columns directly tells us if a 
certain orbit is exceptional or not. In case of coincidence, 
we have left the second column blank. From an experi- 
mental point of view, the appropriate order for the filters 
is something that the experimentalists should elucidate, 
since it is related to the physical substrate used to im- 
plement the qubits, and the resources at their disposal. 
If minimizing the number of two-qubit entangling opera- 
tions is a critical factor, because the fidelity in performing 
such quantum gates is lower than desirable, then the ap- 
propriate order is (\E\,x') : once the number of gates is 
minimized, then it is the turn of reducing the preparation 
depth. 

To complete our results in fl5j we provide, in the two 
final columns for each LC class, an optimum graph result- 
ing from the filters applied in the order min(|i?|, x' ■> #), 
and another one applying the filters as min(x / , \E\, jf) 
(if it is not coincident with the former; in case of co- 
incidence, the second final column is left blank). The 
edges of the graph are listed, with vertices indexed as 
0, . . . , n — 1. Moreover, edges are divided in classes (en- 
closed by parentheses) that define a proper edge-coloring 
(with x' colors). This information is equivalent to pro- 
viding an optimum circuit [with a number of time steps 
equal to x'(G')] f° r preparing the optimum representa- 
tive graph G of the class. The graphs and circuits in 
Fig. [3] have been designed according to the information 
in the supplementary material [Til ] for LC class number 
136. 



IV. ONE-QUBIT GATES 

Assuming that an experimentalist has prepared the op- 
timum representative graph state |G') corresponding to 
the desired state |G), he or she needs to know at least 
one sequence of local complementations connecting \G') 
with \G). The length of this LC sequence is not rele- 
vant, due to the possibility of implementation of these 
local operations as one-qubit gates in only one time step, 
as was already discussed above. However, finding a way 
in the orbit from \G'} to \G) is a hard task. To make 
the entire orbit-based preparation procedure practical, 
we have designed a computer program that accomplishes 
this task. Very briefly, it uses the information about the 
graph G related to the state \G) that one wishes to ob- 
tain. The input is the information about the edges. The 
program finds the optimal graph(s) G' with respect to 
both the number of edges and chromatic index (these two 
quantities are also part of the output), and provides a se- 
quence of LC operations transforming \G') into \G). This 
computer program is included as supplementary material 



V. EXAMPLE 

In order to clarify the whole process, we go back to 
the graph in Fig. Q] As we mentioned, G is a class-1 
graph. Therefore, the preparation depth of \G) using 
only controlled- Z gates is 7. The orbit-based procedure 
allows us to reduce significativcly the preparation depth 
and the number of controlled- Z gates. It consists of the 
following steps: 

(I) To identify the orbit or LC equivalence class 
the graph state \G) belongs to, we calculate the 
cardinality- multiplicity invariants The result is: 



13, 



we 



{0i70, 135, 3i2, 4 7 }. Therefore, after consulting 
conclude that the graph state \G) belongs to the LC class 
number 68 



(II) Also in 15| we find the optimum representative 
graph G': it is the eight- vertex linear cluster LC% (see 
Fig. H]). Graph LCg is a class-1 graph, whose maximum 
degree is A(G") = 2. Hence its preparation depth is 
2, which means a remarkable saving in the preparation 
depth compared to that of \G). 

(III) Therefore, it is worth preparing |G) by prepar- 
ing \G') and then applying suitable one-qubit gates. The 
program in fill ] outputs a sequence of local complementa- 
tions which connects G' to G. For instance, the sequence 
of LC operations applied on vertices 6, 7, 4, 5, and 2 in 
graph G' enables us to obtain G. Denoting the corre- 
sponding scries of local Clifford operations by t(G'), we 
have \G) = t(G')\G'). Applying Eq. (fTJ) and re-arranging 
terms so that r(G') = Ri, where Ri is a specific 
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gate on qubit i, we obtain that 

Ro =cxp(i| ( ri )) , (3a) 

Ri = exp [i^vP^J . (3b) 

i? 2 = exp (-« J4 2) ) exp (^f) , (3c) 

ifc = exp (i|4 3 >) , (3d) 
i? 4 = exp (^i 4) ) exp (-if4 4) ) exp (^ 4) ) , (3e) 

i? 5 = exp (-i J4 5) ) ^P (^af) , (3f) 

i? 6 = exp (i|*< 6) ) exp (-«^4 6) ) > (3g) 

i? 7 = exp exp (-* J<4 ?) ) • (3h) 



In addition, the number of controlled- .Z gates neces- 
sary to get \G) is remarkably reduced (six two-qubit gates 
fewer than in the standard preparation method). The op- 
timum circuit for preparing \G), with a preparation depth 
equal to 3, is the one in Fig. [5] 

VI. CONCLUSIONS 

We have proposed a procedure for the optimal prepa- 
ration of any of the more than 1.65 x 10 11 graph states 
with up to 12 qubits, based on their entanglement prop- 
erties. Optimal means with both (a) a minimum number 
of entangling gates and (b) a minimum number of time 
steps, when possible, or choosing between (a) or (b), in 
the other cases. The preference will depend on the par- 
ticular physical system we are considering. The main 
goal has been to provide in a single package all the tools 
needed to rapidly identify the entanglement class the tar- 
get state belongs to, and then easily find the correspond- 
ing optimal circuit(s) of entangling gates, and finally the 
explicit additional one-qubit gates needed to prepare the 
target, starting with a pure product state and assuming 
only arbitrary one-qubit gates and controlled-Z gates, 
which constitutes the most common scenario for practi- 
cal purposes. 



The resul ts p resented in this paper goes beyond those 
in 0, H, [T3, [l!| : the classification of entanglement for a 
highly relevant family of qubit pure states (graph states 
and, by extension, stabilizer states) of 9, 10, 11, and 12 
qubits. In total, almost 1.3 x 10 6 entanglement classes 
arc introduced. 

The results have experimental relevance. For example, 
imagine an experimentalist in the field of trapped ions 
who wants to prepare graph states. The experimentalist 
knows that he can keep, e.g., nine ions (qubits) isolated 
from external influences for a given period of time, and 
knows that during this time he can perform a maximum 
of m two-qubit entangling operations with an efficiency 
above a certain threshold. The experimentalist wants to 
know which classes of graph states (which classes of en- 
tanglement) are a reasonable target with these resources. 
The results in this paper allow him to answer this ques- 
tion: if m = 8, he can prepare 47 different classes (classes 
147-193 in [lj|); if m = 9, he can prepare 47 + 95 = 142 
different classes (classes 147-288 in [Hi), etc. Moreover, 
the paper tells the experimentalist which is the optimum 
sequence of lasers (gates) required for preparing any state 
of each class. 

More interestingly, consider that the experimentalist 
wants to prepare a specific nine-qubit graph state. The 
paper provides the simplest known protocol to identify 
which entanglement class the target state belongs to, 
and give the simplest circuit to prepare it, where sim- 
plest means in most cases the one requiring the minimum 
number of entangling gates and computational steps or, 
in those cases in which such a circuit does not exist, gives 
a circuit requiring the minimum number of entangling 
gates, and a circuit requiring minimum depth. 
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FIG. 3. Two optimum representatives for LC class number 136 (up). The graph on the left is obtained by applying over the 
entire orbit the filter min(|_E|, x'), i- e -, nrs t minimizing over \E\ and then over \ ■> whereas the one on the right is obtained by 
applying the filter min(x', Any graph state belonging to LC class number 136 could be prepared by means of the circuits 
depicted under those representatives. The circuit on the left prioritizes a lower number of entangling gates; the one on the 
right prioritizes a lower preparation depth. Ri and R[ are one-qubit gates. They are specific for the state of the class 136 we 
want to prepare. 
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FIG. 4. Graph state LCs, optimum representative of orbit 



Qubito — ♦ 
Qubit 1 
Qubit 2 
Qubit 3 
Qubit 4 
Qubit 5 
Qubit 6 
Qubit 7 — ♦ 



FIG. 5. Optimum circuit for preparing the graph state corre- 
sponding to Fig. [T] 



